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We have studied the kinetics of cluster formation for dynamical systems of dimensions up to n = 8 
interacting through elastic collisions or coalescence. These systems could serve as possible models 
for gas kinetics, polymerization and self-assembly. In the case of elastic collisions, we found that 
the cluster size probability distribution undergoes a phase transition at a critical time which can be 
predicted from the average time between collisions. This enables forecasting of rare events based 
on limited statistical sampling of the collision dynamics over short time windows. The analysis 
was extended to L p -normed spaces (p = 1, . . . , oo) to allow for some amount of interpenetration 
or volume exclusion. The results for the elastic collisions are consistent with previously published 
low-dimensional results in that a power law is observed for the empirical cluster size distribution 
at the critical time. We found that the same power law also exists for all dimensions n = 2, . . . , 8, 
2D L p norms, and even for coalescing collisions in 2D. This broad universality in behavior may be 
indicative of a more fundamental process governing the growth of clusters. 
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I. INTRODUCTION 



This paper is a study of the statistical behavior of 
the dynamics of clusters which are allowed to interact 
through elastic collisions or by coalescence. The elastic 
collision dynamics are based on a ballistic billiard model 
analyzed theoretically by Sinai [l|, Q. Cluster growth 
and self-assembly processes are relevant to a variety of 
research fields, including chemistry, materials science, 
physics and earth sciences. The study of such random 
processes can reveal information on the nature of col- 
lective interactions as well as make predictions on the 
occurrence of rare and catastrophic events. Early theo- 
retical studies on cluster dynamics originate in the work 
of Bogoliubov, who showed that in the gas phase, groups 
of particles with short-ranged interactions behave like in- 
dependent clusters @ . Sinai provided a proof of cluster 
dynamics for colliding billiards for one dimensional (ID) 
systems and subsequently for higher dimensions (re- 
stricted to sufficiently low densities). Sinai also proved 
ergodicity of the classical billiard model 0]. The statis- 
tical properties of cluster dynamics has been studied for 
the 2D case with frictionless elastic billiards In this 
paper we extended this statistical analysis of ballistic bil- 
liards to higher dimensions (n) up to n = 8, higher densi- 
ties (p), L p -normed distance metrics (p — 1, . . . , oo) and 
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to the case of coalescing billiards. We have found a high 
degree of universality which suggests that the dynamics 
of clustering are relatively independent of the details. 

In an ensemble of interacting particles we may observe 
a phase transition where a dominant cluster emerges Q ■ 
In a classical Sinai billiard consisting of elastic collisions, 
the phase transition in the empirical density of clusters is 
not necessarily associated with a phase transition of the 
physical system in the traditional sense, such as a tran- 
sition from liquid to solid or gas to liquid as function of 
temperature or pressure. Instead, one observes a change 
in the empirical cluster density - which plays the role of 
the order parameter - as function of time. Thus, it is 
indicative of the dynamics of the motion rather than a 
configurational change resulting from the variation of an 
intensive variable. The collisions between billiards repre- 
sent the interactions between parts of a system, and the 
transition to a dominant cluster that emerges is a man- 
ifestation of the events leading to a major catastrophic 
event. In the case of coalescence, the phase transition in 
the probability density can be associated with a physical 
change in the properties of the system. Recent exam- 
ples of the analysis of phase transitions in probability 
densities include earthquake prediction [5H9I, economic 
modeling [l(| and models of river networks nj. Estab- 
lished premonitory patterns have allowed the modeling of 
events in complex systems to be predicted using observed 
background activity [12]. The prospect of predicting 
or controlling critical events occurring in a dynamic and 
complex environment is of broad interest. 

In the first part of the paper, we expand the study of 
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dynamical phase transitions to higher dimensional elas- 
tic billiards and describe the statistics of the collisions in 
n-dimcnsional L p -normed spaces, with n — 1, ... ,8 and 
p = 1, . . . , oo. The first result which emerges is the exis- 
tence of the phase transition in dimensions greater than 
?i=2, for higher densities, and for different L p -normed 
spaces (p = 1, . . . , oo). In the Euclidean norm case of 
Sinai billiards, the critical time appears to be indepen- 
dent of the dimensionality of the system (n). Instead, 
this critical time solely depends on the average time be- 
tween collisions (r n ). Another notable finding is that the 
empirical cluster distribution at the critical point obeys 
a power law across all dimensions, densities and norms 
with the same exponent. 

In the second part of the paper, we allow the billiards 
to coalesce and form larger clusters. These collisions can 
be analyzed using a binary tree model [Til ] first devel- 
oped to analyze environmental transport in river net- 
works. Coalescence and coagulation are phenomenon 
that are present in many areas of chemistry (T5 , fl4| . 
Theories of coalescence date back to the work of Smolu- 
chowski [H|,[l6{ m the early 20th century, establishing the 
evolution of the concentration, Ck(t), of clusters of mass 
k using a master equation of the form (discrete case): 



dc k (t) 
dt 



= \Y. KiM^y^Y.KjkcM (1) 
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where Kij is the interaction kernel, which is dependent 
on the collision process of i-mers and j-mers. The first 
term predicts an increase in Ck{t) due to coalescence of 
an i-mei and j-mer; the second term deals with the de- 
crease in Cfc(i) due to fc-mers coalescing with clusters of 
different sizes [l3j]. The theory is based on two impor- 
tant assumptions: coalescence upon collision and the ab- 
sence of hydrodynamic interaction between the different 
i-mers [17J . Recent work in the field has led to correc- 
tions to Smoluchowski's equation. Such works include 
film drainage theory [l8| and studies considering a hy- 
drodynamic interaction term (l7l . [l9| . 

The coalescence process we analyze is similar to that 
of Smoluchowski in the sense that there are no inter- 
particle interactions except for coalescence events which 
occur upon contact, and the process begins at t = with 
a monodispcrsc collection of monomers. The model is 
found to exhibit similar properties with regards to uni- 
versality of the phase transition as the elastic model. The 
results suggest that coalescing processes are governed by 
principles similar to that of non-coalescing billiards. 



{A = (x 1 ,X2,...,x n ) : 

< xi < 1; < x 2 < 1; • • . ; < x n < 1} (2) 



The density of billiards within this domain is: 



V n (R)N 



(3) 



where V n (R) is the volume of a n-dimcnsional hyper- 
sphcrc of radius R and Va„ is the volume of the hyper- 
cubic domain. 

Clusters are defined in the elastic collision model using 
the notion of a A-cluster [2(|. Time is measured by the 
variable t, and A specifies an interval of time. A A- 
cluster is a group of billiards which have effected each 
others' kinematics in the previous time interval A. A 
A-neighbor is defined as two billiards at time t which 
have collided during the time interval [t — A, t]. The set 
of all billiards which have interacted and are linked by 
A-neighbor relationships, are called a A-clustcr. A A- 
cluster's mass, M, is the sum of the mass of the billiards 
which make up the cluster. We use the notation M l A (t) 
at time t to represent the mass of the i th largest A-cluster 
(by mass) and NA{t) as the total number of clusters at 
the time t. 

For the coalescing case, clusters are defined as A- 
ncighbors or A-clusters only if they coalesce upon col- 
lision. Each cluster corresponds to one billiard in the 
hypercubic domain. M\ (t) is the mass of the i th largest 
A-cluster at time t and N&(t) is the total number of 
clusters at time t. 



A. Kinetics of Collisions 

We use a ballistic colliding billiard model for both the 
elastic and coalescing cases. We use an a priori method 
to detect collisions in the Euclidean norm case and coa- 
lescing case. When varying the p norm (p ^ 2), we use 
an a posteriori method to detect collisions. Each of these 
arc further described below. 



1. Sinai Billiard 

Elastic hard sphere collisions are enforced for the Sinai 
billiard. Total energy, E, and momentum, mv, of the 
billiards in the system remain constant: 



II. MODEL AND DEFINITIONS 

We start with a Sinai billiard Q in n dimensions, in- 
volving N spheres positioned inside a frictionless hyper- 
cubic domain. The total mass, m, of the billiards in each 
case is 1 . with radius R. The domain is the set of points 
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(4) 



where |vj| 2 = (v\) 2 + ... + {v l n ) 2 . Incident and reflective 
angles arc identical for billiard-wall collisions. 
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where 1 < i < n. Wc then calculate the associated 
unit normal vector rj: 



k—n 
k=l 



V_ 

m 



(7) 



At the time of collision, the initial (i) velocity of 
a is given by v l a = (v l al , . . . ,v l an ). We denote the 
velocity for b analogously. The initial relative ve- 
locity of a and b is calculated and dotted with the 
unit vector, ?), to find the speed, v\ associated with 
the impulse. 



FIG. 



1: Fractional cluster masses M l A as a function of 



time when they were formed for TV = 10 3 , p = 0.001, in 
dimensions as indicated. Each point corresponds to the 
creation of single cluster of fractional cluster mass M l A 

at time t. A phase transition is observed where a 
dominant cluster emerges at a particular critical time 
t c ,n, which is different for each n dimension. Typical 
results are presented for each dimension model, each 
with equal densities and number of billiards. 



In constructing the simulation, an a priori method was 
used to calculate the next time, t next , of collision between 
two billiards or between billiard and wall. The system 
was advanced to this time, t curren t = t next , the velocities 
redefined for colliding billiards, clusters recorded, and the 
next collision time, t ne xt, computed. 



2. IP -Normed Spaces 

We also investigated the effects of using an L p -norm 
(p ^ 2). The motivation for this is to consider particles 
that are of irregular shape and surfaces not unique to 
Euclidean space in collisions. For the application of L p - 
norms in this context, see [2lf. The L p norm is 



l 'r =ViVa- v b) 

3. The impulse, I = (ii, . . . , 7„), is 



(8) 



f „ m a m b i „ 

I = 2 vl-q 

m a + m b 



(9) 



4. Velocities are reassigned for the two billiards to find 
the final velocity (f) after collision. 



vlj =^aj -Ij/m a , 



€- 



= vlj +Ij/m b 



(10) 



In implementing the V norm dynamics, an a posteriori 
method is used where a constant time interval, <5t, was 
used to advance the billiards at each time step. Billiard- 
billiard intersections or billiard-boundary overlaps are 
used to identify collisions. When those collisions are iden- 
tified, velocities are redefined using the method described 
above. The model then proceeds to the next time step. 
The model was verified against the a priori model in the 
L 2 (Euclidean) case in order to find a time step that is 
appropriate for simulation speed. 



\\x\\ p =(\x 1 \P + \x 2 \P + ... + \x n \P) 



p\l/p 



(5) 



When n-dimensional billiards collide in any L p -norm, mo- 
mentum and energy are conserved by reassigning veloci- 
ties as follows: 

1. The vector normal to the sphere's surface at the 
point of collision is calculated between two col- 
liding billiards, a and b. First, by calculating 
the normal component along each dimension, ff = 
(?7i) %)•••) Vn) using the billiard positions ri = 
(r a i, ■ ■ ■ , r an ) and f b = (r b i, r bn ). 



3. Coalescing Billiard 

Coalescing billiards were introduced to model polymer- 
ization reactions and self-assembly processes. In a dif- 



ferent but mathematically similar context, Zaliapin 11 1 
studied transport in river networks. This description of 
coalescing processes can include, for example, emulsions 
of oil in water, or colloidal particles flocculation in per- 
colation analysis [23 | ■ As a motivation for the analysis 
one may consider a growing spherical polymer or colloid 
particle. The probability of coalescence for collisions of 
particles is assumed to follow an Arrhenius law: 



m=\v (r a i - r bi ) p 1 \ 
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(11) 



4 



1 



0) 

N 0.1 r 




rq 1 e-05 r 



c ie-06l ■ — ■ — ■ — 1 

lu 1 10 100 1000 

Cluster Size M 

8D 4D « -2.5 Slope 

6D * 2D □ 

FIG. 2: The power- law empirical cluster size 
distribution at the critical time, t c for TV = 1000 
billiards at a density of p = 10~ 4 , for the different 
dimensions as indicated. One can see the approximate 
power law relationship, with /3 sa 5/2. The magnitude of 
the cluster size, M, is plotted along the horizontal axis; 
the empirical cluster density of cluster size M is given 
on the vertical axis based on averaging from 50 trials. 



where T& is a temperature and E a is the activation energy 
required for coalescence. For individual collisions k is the 
conditional probability of a single coalescence event in- 
volving two billiards in a collision taking place where E a 
is fixed and 71- is the total kinetic energy of the two collid- 
ing billiards. Comparing a uniformly distributed random 
number in the interval {z|0 < z < 1} to k, coalescence of 
the two n-sphcres proceeds if z < k; otherwise, an clastic 
collision occurs. In the latter case, the kinematics of the 
clastic collision are the same as that described for the 
Sinai billiards. In the former case (event of coalescence), 
the two coalescing billiards, designated daughters a and 
b, join to form one billiard, the parent, 7r. In doing so, 
the radius, R n , of the parent is defined in terms of the 
volumes of both initial clusters as to maintain constant 
density of billiards in the domain of the n-dimensional 
system, 

= l(Ra) n + (R b ) n ] 1/n ■ (12) 

Mass is also conserved, with = m a + mj,, where m a 
and mt, are the masses of daughters a and b respectively. 
The resulting magnitude of velocity of the parent billiard 
is defined by conservation of energy: 



conserved. The parent's center, (xi,X2, ■ ■ ■ ,x n ), is de- 
fined as the center of mass of the two daughter billiards. 
If the parent extends beyond the boundary of the system 
upon definition, its center is redefined perpendicular to 
the boundary edge so that it is completely within the con- 
fines of the system. In the event that the parent overlaps 
with another billiard in the system, that billiard under- 
goes a k =1.0 probability collision with the parent in the 
same manner as described above. 

As time evolves and billiards collide and coalesce, we 
obtain a binary tree of coalesced billiards. Each cluster is 
an individual growing billiard on the surface. The total 
number of clusters is designated as Na {t) at time t and 
M\(t) represents the mass of the i th largest cluster. Due 
to the computationally demanding nature of this simula- 
tion, we investigated only the Euclidean distance metric 
(L 2 ) together with the a priori method. 

B. Model Parameters 

Both models were simulated varying N . the number of 
billiards 100 < N < 5 ■ 10 3 ; p, the density 10~ 6 < p < 
10 _1 ; E a , the activation energy satisfies < E a . At time 
t = 0, non-overlapping billiards are randomly placed in 
the volume 



t' = {(xi,x 2 , • • • , x n ) : R < xi < 1 — R; 

...;R<x n <l-R} (14) 

Particles are assigned an initial Maxwellian distribution: 

Several L p norms in the range 1 < p < oo were investi- 
gated. The temperature is held constant at T = 1. 

C. Phase Transition in Cluster Dynamics 

During an interval < t < t c the cluster distribution 
evolves continuously and without gaps. In this regime, 
the size of the largest cluster at any time, (t) , is not 
substantially larger than that of the second largest clus- 
ter, M^(t). A dramatic change (Figure [I} occurs at 
t = t c typical of a phase transition in which a dominant 
cluster appears. The following definition of the critical 
time t c has been proposed Q 



The direction of the parent's velocity vector is defined as 
that which results from a completely inelastic collision 
between the two daughter billiards where momentum is 



t c = M{t : M} > Ml, r > t,i > 1} (16) 

This definition has the operational disadvantage of not 
being a stopping time, meaning that the occurrence of 
the phase transition cannot be decided based on previ- 
ous knowledge of history available until the present time. 
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FIG. 3: Average critical times, t c , plotted for each 
dimension with TV = 1000 varying density as indicated. 



Our results suggest that despite the lack of a stopping 
time definition, perhaps an even more important obser- 
vation is the empirical dependence of t c on the average 
time between collisions, (t„), highlighting a potential way 
to predict its occurrence given some amount of statisti- 
cal sampling. The prediction of t c is important for the 
analysis of earthquake events , interacting billiards 
p3|, ec onomic models [l(| and geographical river mod- 
els [llj. In each of these instances, it has been suggested 
that cluster dynamics could be used to further under- 
stand the system or identify and predict the occurrence 
of a critical event. 



III. RESULTS 

A. Onset of the Phase Transition 

In Fig.[T]we vary the dimension with initial parameters 
p = 10~ 4 , TV = 1000. For each case, n-D, a similar 
characteristic formation of clusters occurs with respect to 
time. The 2D case is in agreement with results previously 
published by Gabrielov et al. For each n a phase 
transition is found at a critical time, t c (n), associated 
with a rapid increase in the growth of a dominant cluster, 
Mi(f c ). 

As the dimension n is varied, while keeping other pa- 
rameters constant, the time scale for the evolution of clus- 
ter dynamics is altered, as seen in Fig. [1] There are two 
contributions leading to this change. 

1. High packing densities become less available for 
packing spheres at higher dimensions. This point 
is discussed by Skoge et al. [24j ]. 

2. The critical time, t c , is found to be dependent only 
on the average time between billiard collisions, (t„) 



for n > 2. This empirical dependence is shown in 
Fig. [3] for approximately equal (r n ) at constant TV. 

The average time between collisions, (r„) is calculated 
from 



(Tn) 



t(M£ 



0.95) 



Z(A) 



(17) 



where Z(A) is the total number of billiard-billiard colli- 
sions that have taken place in the time interval [0, A]. 

For fixed density p the average time between collisions 
(r n ) first increases, then decreases with increasing n, as 
explained by the first point. Due to the critical time's 
dependence on (r„), it scales accordingly. There is a gen- 
eral trend that the dimension, n, for which (r n ) reaches 
its maximum value, increases with decreasing density. 
Thus, we observe similar cluster distributions at each 
dimension, scaled by a factor. From our knowledge of 
the behavior of systems of constant density and varying 
dimension, we have that as the density p decreases, the 
largest critical time is associated with higher dimensions 
n. 

These empirical results provide a predictive tool to es- 
timate the critical time of the system: by estimating 
(r n ) from limited statistical sampling, one may predict 
t c by using a multiplicative factor. This multiplicative 
factor is computed from our data in Fig. @] Data col- 
lected and fit found that for TV = 1000, p = 2, the critical 
time, t c , as a function of average time between collisions, 
< (t„) < 0.45, for 2D data is t c = 0.63 + 471.5(r 2 ), 
for 4D data is t c = 0.24 + 485.22(r 4 ), and 6D data is t c 
= 0.200 + 454.40(t 6 ). This relation for i c ((r„)) must be 
calculated for each value of TV. This gives an expecta- 
tion value for the critical time, i c , expected ((t)) based on 
parameters of the system. The error for actual measure- 
ments from the expectation value increases in magnitude 
proportional to the critical time. The fractional uncer- 
tainty for each data point, e: 



t c , actual — t c , expected 
t c , expected 



(18) 



was calculated. The fractional uncertainties of the points 
form a normal distribution with a standard deviation of 
cr e = 14% for the 2D case, a e = 16% for the 4D case, 
and a e = 15% for the 6D case. Accounting for this error 
in the fit relationships, the values fall within each other's 
standard deviation for points above t — 0.0025. We per- 
formed additional runs with 1.5 < (r„) < 4.0 for 2D, 4D, 
and 6D and found that the previous fits and standard de- 
viations of the fractional uncertainties a e , from the lower 
(t„) data represented well this higher (r„) data. In ad- 
dition, we considered the standard error, comparing the 
mean critical time for sets of ten runs with equal density 
to the calculated standard error in time between colli- 
sions for that set, (t„). As the critical time increases, so 
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FIG. 4: (a) Right: Critical times for each dimension for approximately equal times between billiard collisions, (t„). 
For approximately equal (r n ), critical times are approximately equal in dimensions n > 2. Left: for 2D, 4D, and 6D, 
the critical time, t c , as a function of average time between collisions, (t„); each point represents an individual trial 
for N = 1000 with varying density 0.2 < p < 0.000001, with < (r„) < 0.45. The calculated fit for the 2D data is t c 
= 0.63 + 471.5(t 2 ), 4D data is t c = 0.24 + 485.22(r 4 ), and 6D data is t c = 0.200 + 454.40(r 6 ). This allows us to 
estimate the critical time through limited statistical sampling of (t„). (b)-(d) depict the distribution of errors, e. 

These distributions are approximately gaussian. 



does the standard error of (r„), as we expect from the 
broadening data as the critical time increases (Fig. H}. 

Another method for identifying the phase transition, 
which was proposed in 0, [H| , analyzes the distribution of 
clusters to M^*' as function of time t. In Fig. [2] we 
plot the empirical cluster density distribution, using fifty 
(50) independent trials with identical initial parameters 
and compare the magnitude of each cluster (M^(t); 1 < 
i < Na(£)) at a time t with the empirical probability 
density of cluster size. This empirical density for the 
cluster size distribution evolves over time according to a 
power law with exponential taper Q : 

g t (M) = M-^exp(^-) (19) 



where the critical exponent is (3 s=s 2.5, and j(t) is a time- 
varying function expressing the tendency for exponential 
tail in the sample. As t approaches t c , "f(t) diverges and a 
pure power law emerges. The power-law fit to our results 
is shown by a straight line in Figure [2] The premoni- 
tion of catastrophic events could be done by sampling 
this statistical distribution over time and ascertain the 
approach to a power law. This provides an indicator of 
the imminence of the phase transition. 

Nearly identical empirical cluster distributions are 
found for each dimension n > 2 at the critical time 
(Fig. [2]). As for the 2D case of [||, varying the density 
did not result in changes to this distribution. For each di- 
mension, the cluster distribution evolves according to the 
power law with exponential tail until the critical time, t c 
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FIG. 5: Evolution of cluster distribution for p = 0.001, 
N = 1000 (a) towards critical time and (b) shortly after 
the critical time as indicated for dimensions n = 2,4,6,8 
as depicted in figure. 

(see Fig. [5J. The unique distribution at the critical time 
indicates a critical exponent for the system at the time 
of phase transition. The critical exponent is consistently 
P 2.5 in each dimension studied for elastic collisions. 

The usefulness of this premonitory sign in predicting 
the imminence of a phase transition is determined by the 
speed at which the power law is approached. We do this 
by tracking clusters over time and taking the logarithm 
of both cluster size and density, and approximate the 
form of the cluster distribution (Eq. ITTJ|) as a low order 
polynomial in log(M) at some time t: 

\og(g t (M)) = Kl (f) - /31og(M) - /sa(t) log(A/) 2 (20) 

The time-evolution of the cluster distribution is recorded 
according to Eq. (j2Tj) , fixing (3 = 2.5, and fitting the 
parameters Ki(t),K 2 (t) for each t. A pure power law 
emerges at the critical time, as n 2 (t) — > 0. 

Cluster density distributions gt(M) were created by 
running simulations fifty times for various parameters as 



FIG. 6: Evolution of K 2 (t) averaged for 50 samples, 
each in 4D, N = 1000, p = 0.1. The average critical 
time for the 50 samples was 0.707; the standard 
deviation for the critical times was 0.10. 



reported in Table |U Results for the time-evolution of the 
coefficient n 2 (t) are presented in Fig. [6l The parameter 
is seen to decay approximately exponentially. The 1/e 
point, denoted t*, can be used as a predictive criterion. 
Fitting the k 2 values using n 2 = A\ -cxp[— (t — A 2 )/t*] + 
B\ , we find values for t* tabulated in Table HI The data 
fits in Table U show that a pure power law fit is first ap- 
proached approximately two to four standard deviations, 
at c , prior to the mean, (t c ), for fifty runs. 

TABLE I: Fitting results for the function (21) 
parameters, n is the dimension (an integer), p is the 
density, N is the number of billiards, (i c ) is the average 

critical time, a tc is the standard deviation from the 
critical time, K 2 ((t c )) is the fit parameter at the average 
critical time, t at n 2 (t) w is the time when the pure 
power law occurs. 



n 


P 


N 


(tc) 




Ka((tc)) 


t* 


t at n 2 {t) w 


2 


0.1 


1000 


1.271 


0.1755 


-0.1324 


0.16601 


0.84 


2 


0.01 


1000 


3.1584 


0.3167 


-0.1525 


0.410102 


2.3 


2 


0.001 


1000 


7.8147 


1.0354 


-0.1312 


1.3766 


6.0 


2 


0.01 


2500 


3.258 


0.2491 


-0.153 


0.351295 


2.47 


4 


0.1 


1000 


0.70751 


0.09624 


-0.1339 


0.126561 


0.54 


4 


0.01 


1000 


5.6819 


0.8655 


-0.1287 


0.74666 


4.18 


4 


0.001 


1000 


32.208 


3.931 


-0.1146 


4.77281 


26.6 


4 


0.01 


2500 


7.5752 


0.6641 


-0.1306 


0.680786 


6.175 


6 


0.01 


2500 


1.9977 


0.2280 


-0.125 


0.266022 


1.54 


6 


0.01 


1000 


1.3255 


0.2133 


-0.1147 


0.203752 


1.02 


6 


0.001 


1000 


11.708 


1.516 


-0.1242 


1.86351 


9.6 



Future research may consider the possible existence for 
an upper critical dimension where the /3 = 5/2 is no 
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FIG. 7: Empirical cluster distribution for varying L p 
-norm spaces as indicated for TV — 1000, p = 10~ 1 for 
50 runs. 
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FIG. 8: Empirical cluster density at critical time, t c for 
E a = 0. with densities and dimensions as indicated. 



longer valid. Previous work has found for a particular 
system of coalescing clusters that there is an exact solu- 
tion of the Smoluchowski equation valid for n > n c = 2 
[Hj|. It is possible that the reason for the observed uni- 
versality for n > 2 could be that these models all exist 
above a certain critical dimension and the results we see 
are independent of dimension. 



B. L p Norm Results 

We now look at dynamical transitions in 2D, but us- 
ing the L p norm metric rather than the Euclidean dis- 
tance. There is an approximately equivalent power law 
distribution at the critical time, as shown in Fig. [7] for 
N = 1000 and for varying p, in the range: 1 < p < oo at 
the same density, p. When L p norms are used to check 
for hard sphere overlap, some amount of interpenctration 
of particles is permitted, as can be seen by the fact that 
an L p -norm can exaggerate the importance of dominant 
vector components regardless of the contribution from 
remaining components. Models of granular media based 
on L p norms have been used to study collisions between 
non-spherical particles plj . 



C. The Coalescence Case 

Coalescence dynamics are relevant to many physical 
processes such as nanoparticle and colloidal growth. We 
consider only dynamical systems in which the total en- 
ergy is conserved. Our investigation here focuses on 
whether the trends observed in the Sinai billiard case, 
namely the behavior of the empirical cluster distribution, 
extend to the coalescence billiard case. 

The empirical cluster mass distribution at t c has a 
power law at t c , as seen in Fig. [5] for the model where all 
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FIG. 9: Empirical cluster density at critical time, for 
2-dimcnsional, TV = 1000, p =0.1, at t c for the case of 
varying percentage of effective collisions forming 
coalesced billiard, with varying E a as indicated. 



collisions result in coalescence. This is similar to what we 
observed in the elastic billiards case. For the coalescing 
cluster distributions, we only considered clusters which 
have been involved in at least one coalescing collision (of 
mass M > 2.0/N). Billiards which have not coalesced are 
not considered, as they do not satisfy the A-neighbor re- 
lationship with any other billiards. The coalescing cluster 
distribution plot compares the magnitude of each cluster 
[2 < i < N^s(t)] against its empirical probability. 

As E a increases, t c also increases due to the decreas- 
ing fraction of successful coalescing collisions. For higher 
E a , fewer small clusters are formed and large clusters 
dominate throughout the trial. Figure |9] shows these dis- 
tributions over multiple runs for different values of E a . 
The critical exponent fit is consistent across all models, 
indicating a large degree of universality in the system. 
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One possible reason for this universality is that all mod- 
els are at the critical dimension. Such a behavior has 
been observed by Kangin the case of the Smoluchowski 
equation and n = 2 [131 ] in the elastic case. Our work 
could be further extended to investigate the form of the 
kernel for Smoluchowski equation for this model. 

The analysis of dynamical phase transitions is rele- 
vant to coalescing clusters. Coalescent events have been 
studied recently in nanoparticle growth trajectories . 
This study may be useful in modeling such phenomena. 

IV. CONCLUSION 

We have extended the analysis of dynamical phase 
transitions to higher dimensions, densities, and norms. 



We have also considered two cases of collisions: clastic 
collisions and coalescing billiards, and for the elastic case 
the analysis included the effects of an L p norm. We have 
found universality in the form of a power law for the prob- 
ability distribution of cluster sizes, with the same critical 
exponent describing these systems. In the non-coalescing 
case, the expectation value of the critical time was shown 
to be determined mainly by the average time between 
collisions. These observations can be used to predict the 
onset of criticality and could be used for applications in 
chemistry such as gas dynamics or polymerization. 
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